Massive photons and Lorentz violation 
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responding dispersion relation and tree-level propagator are determined exactly, and some leading- 
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PACS numbers: 11.30.Cp, 11.30.Er, 11.30.Qc, 12.20.-m, 12.60.Cn 



I. INTRODUCTION 

Recent years have witnessed a growing interest in pre- 
cision tests of Lorentz and CPT invariance. This interest 
can be partly attributed to the availability of new obser- 
vational data and the development of ultra-sensitive ex- 
perimental techniques Moreover, minute departures 
from Lorentz and CPT symmetry can be accommodated 
in various theoretical ideas beyond established physics @ 
providing a phenomenological opportunity to search for 
novel effects possibly arising at the Planck scale. 

At presently attainable energies, Lorentz- and CPT- 
violating effects are expected to be governed by effec- 
tive field theory Q. The general framework based on 
this premise is known as the Standard-Model Exten- 
sion (SME) This framework contains the usual 
Standard Model of particle physics and general relativ- 
ity as limiting cases and therefore permits the identifi- 
cation and analysis of essentially all currently feasible 
Lorentz and CPT tests. To date, the SME has been 
employed for phenomenological studies involving cosmic 



and other particle col- 
neutrinos [ll|, precision 



radiation Q, meson factories 
liders @ , resonance cavities [l0|| , 
spectroscopy [l2l |. and gravity 

The SME has also served as the basis for various the- 
oretical investigations of Lorentz and CPT symmetry. 
These investigations have have shed light on the SME's 
mathematical structure [14j , spontaneous Lorentz break- 
down and Nambu-Goldstone modes [l5|, classical lim- 
its of Lorentz- and CPT- violating physics quan- 
tum corrections and renormalizability [13, [l8[ , features of 
non-renormalizable contributions [19], etc. At the same 
time, these analyses have solidified various aspects of the 
SME's theoretical foundation. 

One topic that has remained comparatively unexplored 
concerns Lorentz and CPT violations in massive vec- 
tor particles in the SME: published work [2(| has been 
focused on phenomenological analyses confined to the 
CPT-even Z° and sectors of the minimal SME. How- 
ever, even for this small region in the full SME's parame- 



ter space basic theoretical results, such as expressions for 
the dispersion relation and propagator, are currently still 
lacking. The present investigation reports on progress to- 
wards filling this gap. 

Massive vector particles with Lorentz and CPT break- 
ing are not only of interest for phenomenological stud- 
ies of the heavy gauge bosons Z° and W^. They have 
previously served a valuable tool for investigations of 
the mass-dimension three Lorentz- and CPT-violating 
Maxwell-Chern-Simons term [l8|]. More importantly, 
massive vector fields play a key role for theoretical stud- 
ies involving the photon because they provide a popular 
method for regularizing infrared divergences in perturba- 
tive quantum-field calculations. In fact, this latter appli- 
cation represents the primary focus of this study. But we 
anticipate that our results for the expression of the dis- 
persion relation and the propagator are equally valid for 
the quadratic part of the full SME's Z° and sectors. 

The outline of the present work is as follows. SectionlLTl 
provides the basic ideas behind the construction of the 
model we are studying. A few remarks on the structure of 
the resulting modified Maxwell-Stucckelberg equations 
are contained in Sec. lIIIl Section HVl determines the exact 
dispersion relation and propagator for our model. Sec- 
tion [V] discusses some leading-order results that are ex- 
pected to be useful for practical calculations, and Sec. lVIl 
gives a brief summary of our results. Some supplemen- 
tary material is collected in various appendices. 



II. MODEL BASICS 

Our primary goal is to introduce a photon mass for reg- 
ularizing infrared divergences in perturbative quantum- 
field calculations in a general gauge. This requires a 
smooth behavior of the internal-symmetry structure in 
the massless limit. Thus, the usual Proca term by it- 
self is insufficient, and the mass needs to be introduced 
via, e.g., the Stueckelberg method [2l[. The original ver- 
sion of this method is appropriate for U(l) gauge theo- 
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ries. The method needs modifications for photons em- 
bedded in the U(l)xSU(2) gauge structure of the Stan- 
dard Model [22j], and it generally fails for non-Abelian 
vector fields. The Lorentz-violating generalization of the 
Stueckelberg method, which is discussed in this section 
and the subsequent one, therefore applies solely to the 
SME's QED limit. We do, however, expect the remain- 
ing part of our study, contained in Sees. II VI andlVl to be 
applicable to more general massive vector fields. 

Our starting point is the usual free-photon Lagrangian 
minimally coupled to an external conserved current j^. 
The generalization of this Lagrangian to include arbitrary 
local, coordinate-independent, translation- and gauge- 
invariant physics with Lorentz- and CPT-symmetry 
breakdown can be cast into the following form [6|: 

C 1 = -lF 2 -A.j 

- \F KX {k F )^F^ + ±e KX ^A x (k AF ) K F^ . (1) 

Here, the field strengths and potentials are real- valued 
and obey the conventional relation Fp, v = dy,A v — d^A^, 
and e KX L LV denotes the totally antisymmetric symbol with 
£ oi23 _ Lorentz and CPT breakdown is controlled 
by the quantities k F and k A F, which are given explicitly 
by the following expressions Q 

oo 

n=2 
oo 

{ k AF ) K ^( k AF +1 \ ai --- a2n - 2 d ai ■ ■ ■ d a2 „_ 2 . (2) 

n=l 

Each k F coefficient as well as each k A d F coefficient is 
taken as nondynamical, spacetime constant, and totally 
symmetric in its a indices. The superscript (d) labels the 
mass dimension of the corresponding photon operator, so 
that the unit of the actual coefficient is [GeV] 4 ~ d . 

The next step is to add a mass-type term 5C m for the 
photon to the above Lagrangian (flj . In the conventional 
case, such a contribution is restricted by Lorentz sym- 
metry. In the present situation, this restriction is absent, 
and more freedom in the choice of SC m exists. This ad- 
ditional freedom partly depends on the type of physics to 
be described. For instance, one may wish to model gen- 
eral Lorentz violation for hypothetical massive photons. 
Alternatively, the aim may be to regularize infrared di- 
vergences that often arise in quantum-field calculations 
involving massless photons governed by Lagrangian ((T|). 
In these two examples, the former may allow more addi- 
tional freedom for SC rn than the latter: consider a situa- 
tion in which the coefficients in Lagrangian (fTJ) are such 
that a subgroup of the Lorentz group remains intact. A 
regulator breaking this residual symmetry may be prob- 
lematic, so that violations of the remaining invariant sub- 
group may have to be excluded from 8C m . 

In the present work, the primary purpose for the intro- 
duction of a photon-mass term is to regularize potential 
infrared divergences in quantum-field contexts. In prin- 
ciple, the structure of SC m can then be chosen as sim- 



ple as possible. For example, the conventional Lorentz- 
symmetric Stueckelberg expression would likely suffice. 

However, we introduce an additional set of Lorentz- 
breaking coefficients for reasons outside the present U(l) 
context: certain aspects of the Stueckelberg model, such 
as the dispersion relation and propagator, will turn out 
to be equally valid for the Lorentz-violating Z° and 
bosons. These particles contain not only the equiva- 
lent of the k F and k A p coefficients, but also, e.g., a 
(k^) fJ,l/ mass- type term. For wider applicability in this 
electroweak context, we therefore also include a (k^)^- 
type contribution into our Stueckelberg mass term. With 
these considerations in mind, we implement the Stueck- 
elberg method [2l| by introducing a scalar field </> in the 
following way: 

8C m = \{d^-mA^ v {d^-mA v ), (3) 

where m denotes the photon mass and 

ff" = jf + G» u . (4) 

Here, the G^ v represents the full-SME generalization of 
the minimal-SME's (k^Y v coefficient and is given by: 

oo 

= J2(G^r va i- a ^-*d ai ...d a2n _ 4 , (5) 

n=2 

where each (G^ 2 ™))^" 1 — <*a»-4 j s contracted with an even 
number of derivatives [23j ] . is spacetime constant, sym- 
metric in fi and v, and totally symmetric in a\ . . . «2n-4- 
This definition still contains some Lorentz-symmetric 
pieces at each mass dimension, which can be eliminated 
if necessary. For example, the Lorentz-covariant contri- 
bution ~ rf" contained in {G^)^ u can be removed by 
taking this coefficient as traceless. 

The inclusion of G^ v does not invalidate the Stueckel- 
berg method. At this point, we leave G^ v undetermined. 
We only require it to be small, so that in the limit of 
vanishing Lorentz violation, ff' v approaches rf v without 
a change in signature and rank [24J. In the present U(l) 
context, where 5C m is intended to serve as a regulator, 
specific regions in {k^P, k AF ) parameter space may only 
be compatible with certain definite choices for G^ v , such 
as G^ u = 0, as discussed above. For applications in the 
heavy-boson context, G^ represents an arbitrary physi- 
cal parameter that can only be fixed by observation. 

As in the ordinary case, the resulting Lagrangian C m = 
C-y + SC m changes under a local gauge transformation 

SAfj, = d^x) , Sep = me(i) (6) 

by total-derivative terms. However, in the absence of 
topological obstructions and with the usual boundary 
conditions, the action and thus the physics remain un- 
changed under the transformation (JS]). The next natural 
step then is to select a gauge-fixing condition F[A, (/>] . As 
usual, a multitude of choices for ^[A, <p\ are acceptable. 
We take 

/[A^^^+W, (7) 
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a choice that will turn out to be convenient for our pur- 
poses. Application of the usual Gaussian smearing pro- 
cedure leads to the following gauge-fixing term to be in- 
cluded into the Lagrangian: 



(8) 



where £ is an arbitrary gauge parameter, as usual. The 
associated Faddeev-Popov determinant det(^f) in the 
path integral results in the ghost term 



-F.P. 



-c(d^d v + £m 2 )c, 



(9) 



where c and c are anticommuting scalars. We mention 
that the possibility of introducing additional Lorentz vi- 
olation into the ghost Lagrangian has been studied [25[ . 
We disregard this option in what follows. 

We are now in the position to present our model La- 
grangian C — Cy + SC m + Cg.t, + Cp.p. explicitly: 

C = -\F 2 - A ■ j + \m 2 A^ v A u - ±-(d^A„) 2 

- \<P{d^d v + £m 2 )0 - c{d^ v d v + £m 2 )c 

- \F KX (k F r x ^F^ + \e^A x {k AF ) K F^ . (10) 

It is apparent that the scalar (f> and the ghosts c and c 
are now uncoupled and can be integrated out of the path 
integral yielding an unobservable normalization constant. 
We will therefore disregard these fields in our subsequent 
analysis. 



III. STRUCTURE OF THE FIELD EQUATIONS 

The equations of motion for our Lagrangian (|10[) read 



m 2 f)^ 



A»=f. (11) 



Owing to its underlying antisymmetric structure, the 
term involving F a p in Eq. (|11[) has vanishing divergence. 
Contraction of the field equations with d v therefore re- 
moves the F a p term and places the constraint 



(12) 



on the A^ term, where we have used our earlier assump- 
tion of a conserved source j v . As per definition, f)"' 3 is of 
rank four, so {d a fj a ^ Ap) projects out one of the degrees 
of freedom contained in A^. The ^-dependent Eq. (fT2"|) 
shows that the source j v does not excite (d a r) a P Ap) . This 
degree of freedom therefore is an auxiliary mode. 

Continuing with this decomposition, it is natural to 
define a component 



Af h = A„ 



(13) 



Employing Eq. ([12]). it is apparent that d ll ff lv Af 1 = 
on shell. With Eq. (fT3"]) at hand, we can now substi- 
tute the decomposition of the vector potential A v — 
Af, h - d u {d a f] afi Ap) j \im 2 ) into the field equations (fTT|) . 
Being a gradient, the auxiliary excitation does not con- 
tribute to F a p. By virtue of its equation of motion (|12l) . 
this component also disappears from the A^ term in 
Eq. pip . The zero divergence of Af, h , on the other hand, 
implies that it vanishes when contracted with the £ term 
in Eq. (|lip . Our decomposition then gives 



+ Kf]<=j' (14) 

for the field equations. Note that the auxiliary compo- 
nent has disappeared entirely (F a p only involves Af, h ) 
and that the source f excites the physical degrees of 
freedom in a ^-independent way. 

IV. FEATURES OF THE GENERAL SOLUTION 

In this section, we study general properties of the solu- 
tions of the equation of motion (|lip . This equation holds 
exactly for our Lorentz-violating Stueckelberg photons at 
the classical level. But the model also exhibits numerous 
similarities to heavy gauge bosons with Lorentz violation. 
For example, the Z° and W ± sectors of the SME also 
contain operators of the type k F , k AFl and G; examples 
of these are kw, &2> and k^, respectively. Although the 
linear equation (| 1 1 1) cannot hold exactly for non-Abelian 
gauge bosons, the types of operators quadratic in the 
fields do agree with those for the photon. One can there- 
fore anticipate that Eq. (fTTj) does govern most aspects 
of the tree-level free behavior of Z° and within the 
SME. In particular, the dispersion relation and propaga- 
tor derived below are expected to hold not only for our 
modified Stueckelberg photons, but also for the heavy 
SME gauge bosons at tree level. 

We begin with the plane-wave dispersion relation. To 
this end, we Fourier transform Eq. (Ilip . which yields 

p 2 rf v - j/y - m 2 rf v + Iv^V^PaPfs 



+ 2{k F ) a ^ v p a p - 2i(k AF ) a e a ^p 



Ot{3(IV 



K = -r (15) 



for the equations of motion in p^-momentum space. 
Here, a bar denotes the Fourier transform, and it is un- 
derstood that the replacement d — > —ip has been imple- 
mented in the Lorentz-violating quantities k F , k AF , and 
G. Let us briefly pause at this point to introduce a more 
concise notation that will enable us present many of our 
subsequent results in a more compact form. We define 



= (k F ) a ^p aPp , E» v = (k AF ) a e a ^pp , (16) 

because in the dispersion relation and the propagator, 
k F and k AF will always appear in this form. Moreover, 
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when convenient we abbreviate the contraction of a sym- 
metric tensor with a 4-vector by placing the vector as a 
super- or subscript on the tensor, e.g., G^ = G^ v p v or 
(G 2 )p = G£G" p^p" , etc. We now rewrite the equations 
of motion (fT5j) simply as S{pY v A u (p) = — J M (p), where 

= {p 2 - m 2 W - (l - p»p v - H& v + 2K" U 

-m 2 G^ + ^(G> 1/ + G; P ^ + (17) 

is the expression of the modified Stueckelberg operator 
in our new notation. 

The plane- wave dispersion relation governs source- free 
motion j M = and can therefore be stated as the usual re- 
quirement that the determinant of 5(p)'V vanishes. With 
the results derived in Appendix[3] this translates into the 
equation 

[S] 4 - 6[S} 2 [S 2 } + 3[S 2 } 2 + 8[S}[S 3 } - 6[S 4 } =0. (18) 

Here, 5™ is the nth matrix power of S{p) ll V: and the 
square brackets denote the matrix trace. 

The various trace expressions in Eq. (TT51) can be cast 
into a factorized form: 

£ det(S) = (r u P„Pu ~ im 2 )Q{p) = . (19) 

The {rf^p^Pv — £m 2 ) piece is associated with the aux- 
iliary mode described by Eq. (fl2]) . The £ independent 
factor Q(p) governs the three physical degrees of free- 
dom. Both factors in the dispersion relation (|19p can 
also contain unphysical Ostrogradski-type degrees of free- 
dom (26|, which are introduced because our Lagrangian 
contains higher derivatives [|| . In the presumed underly- 
ing theory, for which the effective field theory (JTUJ) repre- 
sents the low-energy limit, these modes must be absent. 
Consequently, they should also be eliminated from our 
low-energy model (ITU1) . 

An explicit calculation shows that 

Q = (p 2 -m 2 ) 3 + r 2 {p 2 -m 2 ) 2 +r 1 (p 2 ~m 2 ) + r Q , (20) 

where the coefficients rj = rj (p) are momentum- 
dependent coordinate scalars determined by traces of 
combinations of the various Lorentz-violating tensor ex- 
pressions appearing in Eq. (Ti"5j) . They vanish in the limit 
kF,kAF,G — > 0. The explicit expressions for the rj, 
which can be found in Appendix [B] are not particularly 
transparent. Note that in general the physical dispersion 
relation (|20[) does not represent a true cubic equation in 
the variable (p 2 — m 2 ) because the rj(p) are momentum 
dependent. 

The dispersion relation ([!§]) restricts the set of all 
possible Fourier momenta p^ = (io,p) to those associ- 
ated with plane-wave solutions of the free model. Since 
in general hp , k^p, and G contain high powers of p^, 
there can be a corresponding multitude of plane-wave 
frequencies to (p) for any given wave 3- vector p, a fact re- 
flected in the momentum dependence of the We 



remind the reader that most of these are artifacts of our 
cffcctive-Lagrangian approach and must be eliminated. 
Only those wave momenta that represent perturbations 
of the usual Lorentz-symmetric solutions should be inter- 
preted as physical. In any case, a determination of the 
exact roots of the general dispersion relation p^|) appears 
to be unfeasible. However, an exact discrete symmetry 
of the plane- wave solutions is discussed in Appendix [Cj 
the massless limit is studied in Appendix [D] and some 
leading-order results are presented in Sec. [V] 

In the more general case of non- vanishing sources, the 
construction of solutions can be achieved with propagator 
functions. Paralleling the ordinary Lorentz-symmetric 
case, we implicitly define the p^-momentum space prop- 
agator P{jp)^ v via P(p)^ „ S(p) v \ = —ir]\, where we have 
employed the usual quantum-field convention by includ- 
ing a factor of (— i). It is thus evident that the propaga- 
tor is given by P(j>y v = -iS~ l {pY v . With Eq. flXEJ}, 
we obtain an exact, explicit expression for the modified 
Stueckelberg propagator in momentum space: 

-\[S] 2 S + ^[S 2 ]S+[S]S 2 -S 3 ). (21) 

As before, S n denotes the nth matrix power of S(p)^ v , 
and the matrix trace is abbreviated by square brackets. 
We mention that restricting each of the infinite sums in 
Eq. © to their first term and setting both m and G to 
zero yields the limit in which previous propagator expres- 
sions have been considered [27| . 

When the exact tree-level propagator (|2~1T) is Fourier- 
transformed to position space, an integration contour 
must be selected. As in the conventional case, this choice 
depends on the boundary conditions (e.g., retarded, ad- 
vanced, Feynman, etc.). In the present case, a further 
issue arises. As per our earlier assumption in Sec. HH 
the Lorentz-violating terms are to be treated as pertur- 
bations of the usual Lorentz-invariant solutions, and we 
already commented on the need to eliminate spurious 
modes. In the present context concerning the propagator, 
this may, for example, be achieved with a careful choice 
of (counter)clockwise integration contours that encircle 
only the desired poles for each of the two orientations. 
We remark that these issues are absent when the model 
is restricted to terms of mass dimension three and four. 



V. LEADING-ORDER RESULTS 

The absence of compelling observational evidence for 
departures from Lorentz and CPT symmetry in nature 
implies that if the coefficients kF and kAF are nonzero, 
they must be extremely small. The same reasoning holds 
true for G, at least in the context of the Z° and 
bosons [28|. This fact is consistent with the theoreti- 
cal expectation that potential Lorentz and CPT viola- 
tion in nature would be heavily suppressed, for exam- 
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pie by at least one power of the Planck scale. For most 
phenomenological studies and many theoretical investi- 
gations it is therefore justified to drop higher orders in 
fcp, Uafi and G. This section contains a brief discussion 
of a few leading-order results. 

We begin by considering the physical piece Q of the 
dispersion relation given by Eq. (|20p . We remind the 
reader that Q = cannot be considered a true cubic 
equation in the variable (p 2 — m 2 ) due to the p dependence 
of the coefficients ro, rx, and r 2 . We can nevertheless 
employ the expressions Rj, j £ {1, 2, 3}, for the roots of 
a cubic to transform the single equation Q = into an 
equivalent set of three equations: 

p 2 -m 2 = R^uI^J^afM.G). (22) 

Just as for the original expression Q = 0, it seems unfeasi- 
ble to determine the exact dispersion-relation roots from 
the above three equations: the plane-wave frequency 
u)g contained in the 4-momentum p^ — (cop ,p) still 
appears on both sides of Eq. ([22]) . In particular, this set 
of equations (|22[) will in general possess more than six 
solutions, which is consistent with the fact that Q = 
fails to be a true cubic and can still contain a multitude 
of spurious Ostrogradski modes. 

Before continuing, one may ask whether the six physi- 
cal solutions can exhibit undesirable features. A detailed 
analysis of this question would be interesting but lies out- 
side our present scope. However, continuity implies that 
small Lorentz violation leads to small deviations from 
the conventional dispersion-relation branches, so that the 
solutions must in general be well behaved within the va- 
lidity range of the our effective field theory. Exceptions 
from this expectation require non-generic circumstances, 
such as particular parameter combinations. For exam- 
ple, in the minimal SME only a single coefficient — the 
timelike k AF — is known to be problematic under some 
circumstances [|J . At the level of the dispersion relation, 
this is reflected in the presence of complex- valued solu- 
tions. For example, if m = and k AF is the only non-zero 
Lorentz-violating coefficient, we have 

R 2 = -2{kf F ) 2 2\l{kf F Y + (p-kf F y. (23) 

When k AF is timelike, the j = 2 contribution in Eqs. (|2"2")l 
cannot have real solutions wp for small enough p: Suppose 
there were real solutions uip- G K., then the square root 
would be real, and R 2 would be real and negative. For 
the special value p = 0, this would give lu^ 2 = R 2 {p — 
0) < 0, which contradicts our assumption ujp G R. In 
what follows, we disregard such isolated regions in SME 
parameter space. 

The physical solutions have to represent small per- 
turbations to the Lorentz-symmetric dispersion relation 
p 2 — m 2 = 0. Equation (|2"2"|) establishes that the depar- 
tures from the conventional case are controlled by the 
three Rj. In a concordant frame Q, this means that (as 



opposed to the spurious Ostrogradski modes) the phys- 
ical roots are characterized by small Rj that approach 
zero in the limit of vanishing Lorentz breaking. This ob- 
servation suggests solving Eq. (|2"2"|) perturbatively by in- 
troducing a parameter A multiplying the Rj , and writing 
a Taylor expansion for the plane- wave frequencies: 

4* = 4i + Ha + A2 ^2 + • • • , (24) 

where we have identified the label j of the root expression 
with the label a of the plane-wave frequency. Substitut- 
ing (jlHj) on both sides of Eq. (f2"2")l and expanding in pow- 
ers of A yields an expression from which the coefficients 
of the Taylor series ([Ml) can be determined by matching 
the appropriate terms. Solutions to Eq. (|22p are then ob- 
tained by taking A = 1, and the Lorentz-invariant roots 
lj 3 ^ — wig are recovered for A = 0. Thus, this proce- 
dure continuously connects six exact dispersion-relation 
solutions — two for each j in Eq. (|2"2"|) — to the Lorentz- 
symmetric mass shell. This behavior is consistent with 
that of the physical roots, so that the procedure simulta- 
neously serves as a filter for rejecting the spurious modes, 
which cannot remain finite in the A — > limit. 

A useful practical way to solve for the coefficients in 
Expansion (|24[) is iteration: substitute as a first iteration 
up^ — on the right-hand side of Eq. (|2"2"]l . which 

determines the improved value cji^ for the plane- wave 
frequency on the left-hand side of this equation. Repeat- 
ing this process yields 

<45+i = + m2 + ^'H&P> ^f, k F , G) , (25) 

where we have denoted the iterative step by the sub- 
script n on wi^ 1 . One can check that the nth iterative 
step possesses the correct value for (at least) the first n 
coefficients in Eq. (f24|) . In this process, there is actually 
no need to introduce explicitly the parameter A, as we 
can take it equal to unity immediately. 

Note that particular care may be required for spe- 
cial regions in momentum space. One example concerns 
very small m and p. The exact branches associated to 
the conventional positive- and negative- valued mass shell 
could then lie closely together, so that convergence to 
the correct branch must be ensured. Moreover, the ex- 
act solution for a positive branch may briefly dip into 
the negative spectrum and vice versa leading to an un- 
equal number of positive and negative roots in this small 
momentum-space region. The signs of the square root 
in the recursive equation (|25[) may then have to be ad- 
justed correspondingly. Other issues, such as isolated 
momentum-space points at which the A expansion of s\ 
in the Eq. (|B7[) starts at 0(A), may be resolved by em- 
ploying a different perturbation scheme. 

Since j G {1,2,3} and since for large enough p there is 
a positive and a negative square root for each j, there are 
six equations corresponding to six seemingly independent 
physical plane- wave frequencies for a given p. However, 
the behavior of the modified Stueckelberg operator S(p) 
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under the replacement p — > —p leads to a correspondence 
between the positive- and negative-root solutions. This 
implies that there are in fact only three independent po- 
larization modes, as expected for a massive spin-1 field. 
More details regarding this line of reasoning can be found 
in Appendix [Cl 

Next, we consider the momentum-space propagator 
given by Eq. (|21[) . This propagator describes two phys- 
ical features. First, the poles of the denominator select 
the plane-wave momenta. In a Feynman-diagram con- 
text, for instance, Lorentz- violating corrections to the 
usual poles can cause the propagator to go on-shell, pos- 



sibly allowing processes that are forbidden in the Lorentz- 
symmetric case. Second, the numerator governs features 
associated with the polarization states. For example, 
Lorentz-violating numerator corrections could yield sup- 
pressed processes that would be forbidden by angular- 
momentum conservation in the corresponding Lorentz- 
invariant situation. 

If corrections to both of these features are to be re- 
tained, the denominator and the numerator in Eq. (|21[) 
should be expanded separately to leading order. With 
such an approximation, the dominant contribution to the 
propagator can be separated into the two pieces: 



u 2 2s (f~rrf_ 
-i(p — m ) 



2[K] + G p p - m 2 [G])ri^ - 2K» U + 2i& v + m 2 G^ + p p G v p +p u G^ 



n,=i [p 2 -m 2 - RjWP + m 2 ,p, k AFl k F , G)] 



(P 2 



1 )(P 2 



m 



2[K] -m 2 [G]) 



{fj^PuPv - £ m2 ) IL=i [p 2 -m 2 - R]Wp 2 + m 2 ,p, k AF , k F , G)] 



P»P V 



higher order. (26) 



In the denominators, we have used the leading-order ex- 
pressions for the Rj implementing our previous results 
for the dispersion relation. 

We note that the first piece of the propagator expres- 
sion is independent of £ and only exhibits poles corre- 
sponding to the three physical modes. The second piece 
is £ dependent and contains the additional pole associ- 
ated to the auxiliary mode. As expected, the second piece 
is without physical effects: contraction with a source j u 
yields an overall factor p ■ j, which vanishes for conserved 
currents. In any case, all contributions from the sec- 
ond piece would necessarily have to be proportional to 
(p-j)p^i which is pure gauge. Moreover, the second term 
vanishes identically in Feynman-'t Hooft gauge £ = 1. 
We further remark that the expression (|26l) does not 
propagate the additional multitude of unphysical higher- 
derivative modes discussed in the previous section. This 
is required for a smooth behavior in the limit of van- 
ishing Lorentz breakdown. Some remarks regarding the 
massless limit of the propagator (f2l)|) arc contained in 
Appendix [D] 

In many situations, Lorentz-violating corrections to 
the pole structure of the propagator can be disregarded. 
Then, the expression for propagator can be simplified 
even further with an approximation that essentially cor- 
responds to treating the Lorentz-violating pieces in the 
Lagrangian as interaction terms. This alternative ap- 
proximation is determined next. 

We begin by considering the Lorentz-invariant part Sq 
of Stueckelberg operator 

So(pT v = (P 2 - m 2 ) rT - f 1 - |) PV ■ (27) 
The expression for the corresponding momentum-space 



propagator is given by [2 



rf v — p^p v lm 2 p^p u lm 2 
i t; ^ i ■ 



(28) 



p" — m* p" — £m 2 

Starting from these zeroth-order expressions, we may 
now decompose the full Stueckelberg operator as 



(29) 



with 



8S{pY" = 2i£»" - 2K"" + m 2 G^ v 



--(g^ p v + g;p" + g^g;). (30) 

Suppressing Lorentz indices and the dependence on p for 
brevity, we can now expand the full propagator P about 
its Lorentz-symmetric value Pq as follows: 

P = -i(S - SS)- 1 

= -i[(l-6S So"" 1 ) So] -1 

= -i5,7 1 (1-SS ST 1 )" 1 



P (l + J2(iSSP ) n 



= P + P i5S P a + P iSS P iSS P 



(31) 



where we used Pq 



-iSn 



The infinite geometric se- 



ries should be convergent for small enough Lorentz vio- 
lation SS. In the context of quantum-field perturbation 
theory, the expansion pip can be represented diagram- 
matically with a propagator insertion iSS into internal 
photon lines, as shown in Fig. Q] This result appears 
to be consistent with Matthews' theorem 30]. Like the 
previous expression (|2"6"|) . the expansion (I5T1) of the prop- 
agator around its Lorentz-symmetric value automatically 
accomplishes the elimination of the unphysical poles. 
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FIG. 1. Lorentz-violating propagator insertion. The wavy 
lines represent the usual Lorentz-symmetric Stueckelberg 
propagator given by Eq. ()28|) . The square box denotes the 
the Lorentz-breaking insertion determined by Eq. (|30[l . 



VI. SUMMARY 



In this work, we have employed the Stueckelberg 
method to introduce a mass term into the full SME's 
photon sector in an R^-type gauge, such that the usual 
smooth behavior of the model in the m — > limit is 
maintained. We have discussed the resulting equations of 
motion in the presence of external conserved sources, and 
studied their solutions. Our results include the exact dis- 
persion relation and propagator for free massive photons 
incorporating all possible Lorentz-violating, translation- 
and gauge-invariant Lagrangian contributions at arbi- 
trary mass dimension. From these exact results, we have 
extracted leading-order expressions for the dispersion- 
relation roots and the propagator. The method for ob- 
taining the roots is recursive; it also yields higher-order 
corrections, and it naturally separates the physical solu- 
tions from spurious modes arising from higher-derivative 
terms in the Lagrangian. Finally, we have verified that 
the number of physical modes is three, as expected for a 
massive vector particle. 

Our results arc intended primarily to provide a flexi- 
ble tool for regularizing infrared divergences in quantum 
electrodynamics in the presence of general Lorentz viola- 
tion. A follow-up work in this context is already in prepa- 
ration 31]. In addition, our expressions for the dispersion 
relation and propagator also hold for the quadratic part 
of the SME's heavy-boson sector, and are therefore ex- 
pected to find applications in theoretical and phenomeno- 
logical studies of electroweak physics with Lorentz break- 
ing. In the to — > limit, our results yield the previously 
unknown massless photon propagator in the full SME, 
a quantity indispensable for many Lorentz-symmetry in- 
vestigations in electrodynamics. 
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Appendix A: Matrix relations 

This appendix provides a collection of various textbook 
linear-algebra results necessary for the line of reasoning 
in the main text. For completeness, we have sketched the 
derivation of these results. 

Consider an arbitrary 4x4 matrix A and denote its 
four (not necessarily distinct) eigenvalues by Xj, where 
j = 1, ... ,4. The characteristic polynomial C(A) of A 
can then be expressed as 



C(A) = (A - Ai)(A - A 2 )(A - A 3 )(A - A 4 ) 



(Al) 



By the Cayley-Hamilton theorem, the matrix A must sat- 
isfy C(A) = 0. Expanding the right-hand side of Eq. (|Aip 
and replacing A — !• A then yields 



(A2) 



A 4 + c 3 A 3 + c 2 A 2 + ciA + c l = , 
where the coefficients Cj are given by 
Co = Ai A2A3A4 , 

ci = — A1A2A3 — A1A4A3 — A2A4A3 — Ai A2A4 , 

C2 = Ai A2 + A3A2 + A4A 2 + Ai A3 + A1A4 + AsA 4 , 

c 3 = -Ai-A 2 -A3-A 4 . (A3) 

The next task is to express the Cj in terms of man- 
ifestly component-independent expressions involving A. 
To this end, we transform A to its Jordan normal form 
A' = M~ 1 AM, where M is an appropriate 4x4 matrix. 
Since such a similarity transformation leaves unchanged 
the eigenvalues, it is straightforward to verify that 

det(A') = AiA 2 A 3 A4, 

Tr(A'™) = A? + A2 1 + A 3 l + X2 , (A4) 
where 77 6 N. Inspection now establishes that 
c = dct(A') , 

ci = ±Tr(A ,2 )Tr(A') - ±Tr(A' 3 ) - ±Tr(A') 3 , 
I [Tr(A') 2 - Tr(A' 2 )] , 



(>2 - 
C3 : 



-Tr(A') 



(A5) 



Since both the determinant as well as the trace of a 
square matrix remain unchanged under similarity trans- 
formations, the relations in Eqs. (| A4|) and (|A5j) hold, in 
fact, also for our original 4x4 matrix A. Employing these 
results in Eq. (|A2|) gives 

A 4 -[A] A 3 + ±([A] 2 -[A 2 ])A 2 

+ (i[A 2 ][A] - i[A 3 ] - ±[A] 3 ) A + det(A)l = . (A6) 

Here, we have abbreviated the matrix trace by square 
brackets. This result establishes that the set of matrices 
{1, A, A 2 , A 3 , A 4 } is linearly dependent. 

Equation (|A6j) can be employed to determine further 
relations involving the 4x4 matrix A. For example, 
taking the trace of Eq. (IA6|) yields the expression 

det(A) = ^[A] 4 -i[A 4 ] + i[A 2 ] 2 

+ i[A][A 3 ]- i[A] 2 [A 2 ] (A7) 
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for the determinant of A. Another relation can be ob- 
tained by multiplying Eq. (|A6|) with the inverse matrix 
A -1 . The resulting equation can then be solved for A -1 : 



1 



det(A) 



(i[A 3 ]l + i[A] 3 l-i[A 2 ][A]l 



[A] 2 A + ±[A 2 ]A + [A] A 2 - A 3 



(A8) 



In Eqs. (| A7|) and (IA8|) we have again abbreviated the 
matrix trace by square brackets. 



Appendix B: Dispersion relation 

To present ro, ri, and r 2 appearing in the physical 
dispersion relation (|20j) . we introduce the following more 
efHcient notation. We suppress the AF subscript on Uaf 
and drop the caret: 



(k AF y 



(Bl) 



It is also convenient to extract the traceless part G^ v 
from the full G pv : 



G_p v = G^ v - \ [G\ff v . 



(B2) 



Here, we have again denoted the matrix trace by square 
brackets [G] = G^. Inspection of our model La- 
grangian (fT0|) reveals that for the photon the decompo- 
sition (IB2|) can be implemented by the replacements 



m ->• rh = m(l + |[G] 



my 



(B3) 



where it is understood that £ -1 / 2 has to be placed at the 
appropriate position in the gauge-fixing term. We also 
remind the reader that we denote the contraction of a 
4- vector with a symmetric 2-tensor by writing the vector 
as a superscript or a subscript on the tensor: for example 
(G 2 ) k = G$G% Pa (k AF )*>, etc. 

With this notation, the r coefficients in the physical 
dispersion relation (|20p take the following form: 



l[K} 3 + ^[K 3 ] - 4[K}[K 2 } + 8m 2 K k + 2[K} 2 6p - 2[k 2 ]CF p - 4(fc • p) 2 G£ ~ Am 2 [GK 2 } + 4m 2 [K}[GK] 

- 2m 2 (GKG) p p + 2m 2 [K}(G 2 )P + 2m 2 [GK]GP + 8m 2 (k ■ p)G k + 2m 4 [G 2 K} - m 4 [k][G 2 } - 4m 4 G k k 
+ rh 4 (G 3 ) p p - ±m 4 [G 2 ]G£ - ±m 6 [G 3 ] - 8[Gi^ 3 ] + 8[k}[GK 2 } - 4{k} 2 [GK} + 4[k 2 }[GK} + 8K%6* 

+ Arh 2 [G 2 K 2 ] - 2to 2 [G^] 2 + 2m 2 [(GK) 2 ] + m 2 [K] 2 [G 2 ] - m 2 [K 2 ][G 2 ] - Am 2 [k][G 2 K] + 4m 2 (G k ) 2 

- 4to 2 G^G£ - 2m 4 [G 3 ^] + m 4 [GK}[G 2 } + § fn 4 [k][G 3 } - |rn 6 [G 2 ] 2 + {m 6 [G 4 ] , 



ri = 2[^] 2 -2[K 2 } -A(k-p) 2 
r 2 = 2[k} + G p p + 4k 2 . 



■ 4m 2 k 2 + 2[K]GP + 2m 2 [GK] + m 2 {G 2 ) p p - \m 4 [G 2 ] + 8K% + 4k 2 G p - 4m 2 G k k 



(B4) 



The above dispersion relation can be shown to reduce in the appropriate limit to mass-dimension three results quoted 
in Ref. [3]. 

To determine the quantities Rj appearing in Eq. (|22[) , we apply the usual Tschirnhaus-type transformation for cubic 
equations to Expression (|2"D|) . This yields the equivalent depressed form of the physical dispersion relation Q = 0: 



I 2 ~ 2 

[p — m 



[k] + \G}> + Ik 2 



, 2 -2 

S\ I p — m 



4;„2 



so = , 



where coefficients sq and S\ are given by 

So = # [k] 3 + f [k 3 ] - f [k] [k 2 ] + i(k- P ) 2 [k} + 8m 2 k k k 



8^,2. 2 



k 2 [K] + l[K] 2 G p - l[K 2 ]G p - l\K]{GF p ) 2 

- §(fc • pfQP + ^(Gp 3 - 4m 2 [GK 2 ] + lm 2 [k][GK] - 2m 2 {GKG) p p + §m 2 [£](G 2 )£ + ^m 2 [Gk]& p 

- im 2 GP(G 2 )P - fm 2 fc 2 G£ + 8m 2 (k-p)G k + 2m 4 [G 2 K] - |m 4 [X][G 2 ] - Am 4 G k k - ±m 4 [G 2 ]G£ 
+ m 4 (G 3 )^ - ±m 6 [G 3 ] + ffc 2 ^] 2 + \k 2 [K 2 ] - f[K]K k + fk 2 (k-p) 2 ~ lk 2 (G p ) 2 + 8[K][GK 2 ] 

- 8[GK 3 ] - A[k] 2 [GK] + A[k 2 ][GK] + fk^G p - fk 2 [k}G p - 4m 2 [X][G 2 i^] + 2m 2 [(dftT) 2 ] 
+ 4m 2 [G 2 i^ 2 ] - 2m 2 [GK} 2 + m 2 [K} 2 [G 2 } - rn 2 [K 2 ][G 2 ] - §m 2 fc 2 [G£] + 4m 2 (G^) 2 



fm 2 k 4 + lm 2 [k}G_l - lm 2 k 2 (G 2 )P p - 2m 4 [G i K] 



V [Gk][G 2 ] + lm 4 [k][G 



- lm 2 G k k GP 
\m 4 k 2 [G 2 } 



|m 6 [G 2 ] 2 + im 6 [G 4 ] + f k 4 [K] - fk 2 K k - fk 4 G? p + fm 2 k 2 G k + ^k 6 . 



Sl = i[K] 2 -2[K 2 ]-A(k-p) 



4 1- ) 9 

2 1 4m 2 fc 2 



[K]G_p - i(Gp 2 + 2m 2 [GK] + m 2 (G 2 )P p - ^m 4 [G 2 ] + 8K k 



^k 2 [K} + U 2 G_P-Am 2 G k k 



16 7,4 
— K . 



(B5) 



(B6) 
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With So and Si at hand, the usual formulae for the roots of a cubic determine explicit expressions for the Rj 
quantities. These expressions can then be used for the iteration Eq. ([251) : Suppose we want to calculate the first-order 
frequencies from the zeroth-order solution {\/p 2 + m 2 ,p), which corresponds to the first iteration step in Eq. (|25l) . 
With (\/p 2 + m 2 ,p) and Eqs. (|B6|) . we can in principle check the sign of the discriminant 4si 3 + 27s 2 , at least 



to leading order. A positive discriminant is associated with complex- valued Rj, which would appear to lead to 
unphysical roots. We therefore take the discriminant to be non-positive in what follows. The Rj can then be 
expressed trigonometrically: 



Ri =2\ cos 



- cos 
3 



3so 
2si 



3 \ 2tt 



-§[x]-§<2£ 



|fc 2 



je {1,2,3} 



(B7) 



These explicit Rj can now, for example, be employed in 
the context of Eq. ([25]) . 

We remark that for the special case of a vanishing dis- 
criminant, Eq. (jB7[) determines a multiple root, so that 
at least two physical excitations propagate without bire- 
fringence. We also note that for non-positive discrim- 
inant the argument of the cosine function in Eq. (|BT|) 
is real, and its magnitude therefore remains bounded by 
one. The limit of zero Lorentz violation is then deter- 
mined by st, [K], G_p, and k 2 . These coefficients will 
generically be momentum dependent. In particular, the 
zero components for the unphysical solutions may diverge 
in this limit, so that the corresponding behavior of the 
Rj is unclear a priori. However, for the physical solu- 
tions the plane-wave frequencies must remain bounded 
because they are defined to approach the conventional 
expressions in the Lorentz-symmetric limit. For vanish- 
ing Lorentz breaking, the physical Rj are therefore seen 
to approach zero, as required by self-consistency. 



Appendix C: Discrete symmetries 

An interesting question concerns the number of inde- 
pendent plane-wave solutions. Some insight into this is- 
sue can be gained by investigating discrete symmetries 
of both the plane-wave polarization vectors and the dis- 
persion relation fp~9]) . 

Since the Lagrangian (jlOj) is real, the momentum-space 
Stueckelberg operator must be hermitian 

S""(p)* = S^(p), (CI) 

which can also be established directly by inspection of 
Eq. ([T7| . Next, we consider S^ v (p) under p — > —p. Only 
the k AF piece reverses its sign because it is odd in p. All 
other terms are p even and remain unchanged. We also 
need the behavior of S tll/ (p) under complex conjugation. 
Since the position-space Stueckelberg operator is purely 
real, complex-valued contributions to the momentum- 
space S^" (p) can only arise from the Fourier replacement 
d — > —ip. We will focus on real-valued plane-wave mo- 
menta p since otherwise an interpretation as a propagat- 
ing solution is difficult. It then follows that the p-odd 
k AF term is purely imaginary, while all other terms are 



even in p even and thus real. Complex conjugation of 
S ilv {p) therefore only changes the sign of its k AF piece. It 
thus becomes evident that the momentum-space Stueck- 
elberg operator is left unaffected under the combination 
of the two operations: 

S^i+p, +k AF , k F , G) = S^i-p, +k AF , k F , G)* 

= S^(- P ,+k AF ,k F ,G), (C2) 

where the Hermiticity of S^ v has been used. Since any 
determinant is invariant under transposition, we find 

det S(+p, k AF , k F ,G) = det S(—p, k AF , k F , G) , (C3) 

i.e., p — > —p represents a discrete symmetry of the dis- 
persion relation. This invariance also becomes clear when 
Eqs. (Jinj), (HU, and CE} are combined. 

Continuing our analysis in a concordant frame we 
fix a wave 3-vector p and imagine solving the disper- 
sion relation ([TO]) . Following the reasoning presented in 
Sec. El we can eliminate the spurious Ostrogradski modes 
and focus on the six physical branches of the dispersion- 
relation solution with exact uP^ . In these expressions, 
the ± label must correspond to the sign of the solution, 
a property that follows immediately from the required 
proximity to the Lorentz-symmetric sheets 32] . Without 
loss of generality, we now consider the specific plane- wave 
4- momentum with positive energy (p 3 + )^ = (^°l^,—p)- 
The p — > —p symmetry of the dispersion relation im- 
plies that (p J _)^ = (—oplp,p) also represents a solution, 
but one with negative frequency. Moreover, this solution 
must be physical: since p J + is close to the conventional 
positive sheet, p 3 _ must be a perturbation of the negative 
sheet. With a suitable selection of the j labels, we may 
therefore conclude uj j ~~~ — —uj^i. In summary, the set of 
physical plane- wave momenta is given by 

{(p!)" = (<4' + >£> . W-Y = . (C4) 

and is thus determined entirely by the expressions for the 
three positive- frequency solutions 

In addition to the plane- wave momenta (pj_) M , the as- 
sociated polarization vectors {e±Y(p) are needed to con- 
struct a plane wave. Given the £ 3 + {p) associated with the 
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positive-frequency roots, we can construct the remain- 
ing negative-frequency s 3 _{p) as follows. Up to normal- 
ization, the polarization vectors are by definition in the 
kernel of the Stueckelberg operator: S , p) e J ±(p) = 0. 

This definition contains S^JZ, —p) e J + (—p) = as a spe- 



cial case. Complex conjugation combined with Eq. (|C2j) 
then yields S(— uj J ^t,p) £+*(— p) = 0. Together with the 
plane- wave momentum property (|C4|) . this equation is 
identified as the defining relation for e 3 _(p). Using suit- 
able normalizations, we arrive at 



(C5) 



We remark that this relation continues to hold for 
dispersion-relation roots with higher multiplicities when 
the bases of KerS'(— are chosen appropriately. 

An arbitrary, free, real- valued, physical solution A^(x) 
can be represented by the following Fourier superposition 
of plane waves pjjjj : 



e -ip J + -x 



d 3 p 

+ al(p) {e 3 _Y(p) e - ip -- x ] + c.c. (C6) 



This expression suggests that six Fourier coefficients 
{a° + {p) 1 a J _{p)} with j = 1,2,3 must be specified inde- 
pendently to determine a general solution. However, a 
change of integration variable p — > —p in the expression 
associated with the negative-frequency solutions together 
with Eqs. ([Cll and (|ClI) yields 



A»(x) 



d 3 p 



g hi, -:\ : - I 



¥ (f)(A_r (p) e-^ + cc, (C7) 



where we have defined b J (p) = a 3 + (p) 



~3 * 



(-p). This 



establishes that the six physical dispersion-relation roots 
are associated with three degrees of freedom. 



Appendix D: Massless limit 

In the context of regularizing infrared divergences in 
quantum-field perturbation theory, the massless limit of 
our Lorentz-violating Stueckelberg model is of particu- 
lar importance. The photon mass was introduced into 
our Lagrangian via the Stueckelberg method, which en- 
sures the proper behavior in this limit. Nevertheless, it 
is interesting to study explicitly the dispersion relation 
and the propagator of our model for vanishing m, for 
example, to compare with previous results for massless 
Lorentz-violating photons. 

We begin by defining the contraction 

J = e aMS e KX ^p pPa p T p u {k F ) a ^{k F ) x ^{k F y s ^ (Dl) 



of the Lorentz-violating kp. One can show that the co- 
ordinate scalar I obeys the following identities: 

T2 IG% = 2[GX 3 ] - mK 2 ][K] + [GK]([K] 2 - [K 2 ]) 



\i P 2 = z[k][k 2 ]-[kf ~2[k*]. 

With these identities, the full massless dispersion rela- 
tion (IT51) can be cast into the form 



(D2) 



i&et S = (ffp^) 2 (p 4 + 2[K]p 2 

- 4(fc • p) 2 + Ak 2 p 2 + 2[K} 2 - 2[K 2 ] ) . (D3) 



8K k k 



It thus becomes apparent that in the massless limit the 
expression (t)^ v ) can be factored out of Q. 

The next step is to identify the physical part of the 
massless dispersion relation (ID3|) . Clearly, one of the 
{fj^PfiPv) factors originates from the auxiliary mode and 
remains unphysical. To gain insight into the other fac- 
tors, it is instructive to recall the model Lagrangian (|10p . 
In the zero-mass limit, G enters the equations of motion 
for the photon only via the gauge-fixing term £ g .f. and 
should therefore not lead to observable effects. It then 
follows that the second {fj^p^Pu) factor must also be un- 
physical because it contains fj = rj + G. We conclude that 
the last factor in Eq. (ID3[) . which is free of G, governs the 
physical excitations for m = 0. This factor is in agree- 
ment with the general photon dispersion relation given 
in Ref. [1], as required by consistency. 

The exact massive propagator (|2"Tj) as well as its 
leading-order approximation (|26p exhibit six physical 
poles that must be encircled by the integration contours. 
The above argument identifies those two of these six poles 
that should become unphysical in the massless case. The 
question then arises as to whether the unphysical char- 
acter of these two modes is reflected in the pole structure 
of the propagator in the m — >• limit. This may be es- 
tablished by demonstrating that the (jf v PiiPv) factor in 
the denominator of the propagator is cancelled by a corre- 
sponding (ff ,v p uPv) factor in the numerator for vanishing 
mass. An investigation of this cancellation for the exact 
expression (|2"Tj) would be interesting but lies outside the 
scope of the present work. We instead focus on the on 
the leading-order piece (|2"B1) of the propagator, which is 
appropriate for the majority of practical applications. 

In the discussion of the leading-order propagator (|26|) , 
it was already mentioned that the (1 — £) contribution in 
the second line leaves unaffected observable quantities. 
We may therefore focus on the physical contribution con- 
tained in the first line of Eq. (|2^|) . In the m — > limit, 
the numerator —iN^ v of this first piece has the form 

= p 2 {2\K\^ v - 2K» V + 2i£f" / + p^G v v + p v G%) 

+ P 2 {p 2 +G^. (D4) 

The rf tlJ term in the second line of the above Eq. (|D4|) 
contains the factor p 2 + G v p = rf v p il p ll . All terms in the 
first line represent leading-order corrections, which re- 
main unaffected by the addition of higher-order Lorentz- 
violating terms. A judicious choice for such higher-order 
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contributions is p 2 —> p 2 + — fj^p^pu in the overall p 2 
factor. It thus becomes evident that in the physical piece 
of the leading-order propagator, the factor {ff u p^Pv) ap- 
pears in both the numerator and the denominator and 
therefore cancels. We conclude that the poles that be- 
come unphysical in the massless limit disappear from the 
propagator, at least to leading order. 

Since G arises in the massless case solely from the 
gauge-fixing term, measurable quantities are necessarily 
free of G, even though p^&i and p^Git terms still ap- 



pear in the numerator of the propagator (I26[> and, equiv- 
alently, in the propagator insertion (|30|) . When the prop- 
agator is contracted with a conserved current ] u , G in- 
dependence can be seen explicitly as follows. For the 
p^Gp term, such a contraction yields an expression pro- 
portional to which is pure gauge and can therefore 
not lead to observable effects. For the p v G^ term, the 
contraction produces an overall p ■ j, which vanishes due 
to current conservation. It is thus clear that G is indeed 
unobservable in this context, as expected. 
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